Appearance
Functional genomics analysis reveals the evolutionary adaptation and demographic history of pygmy lorises
Summary
Problem: 链猴亚目（strepsirrhines）在灵长类系统发育中占据关键节点，但高质量基因组资源稀缺，阻碍了对该亚目适应性进化机制的研究。慢懒猴具有低代谢率、缓慢运动、冬眠、毒素抗性等独特性状，但其遗传基础此前基本未知。 Approach: 作者组装了倭蜂猴（Nycticebus pygmaeus）的染色体水平基因组，并对50只倭蜂猴和6只孟加拉蜂猴（N. bengalensis）进行全基因组重测序。整合比较基因组学（基因家族扩张/收缩、正选择、趋同进化分析）、转录组分析、功能实验（酶活性测定、细胞转染、免疫共沉淀）以及群体基因组学（PSMC/MSMC）方法，系统探究适应性进化的遗传机制和种群历史。 Finding: ①慢懒猴中多个解毒相关基因家族特异性扩张，其中GSTA基因家族在肝脏中高表达且拷贝数显著多于其他哺乳动物；②倭蜂猴与考拉在PITRM1基因上存在趋同氨基酸替换，导致酶活性降低约两倍，可能与低代谢率相关；③MYOF基因受正选择，功能实验显示其上调快肌纤维基因（MyHC2a、MyHC2x）的能力减弱，与慢懒猴慢肌纤维占优的表型一致；④PER2基因受正选择，S1271A突变削弱了PER2与CRY蛋白的结合能力，可能与异常昼夜节律有关；⑤群体基因组分析显示倭蜂猴与孟加拉蜂猴在约100万年前分化，此后有效种群大小呈相反变化趋势，暗示种间竞争，且两者在过去6万年间种群规模大幅下降。 Significance: 该研究填补了链猴亚目高质量基因组资源的空白，首次从分子水平解释了慢懒猴多种独特适应性性状的遗传基础，为理解灵长类适应性进化提供了新视角，同时为慢懒猴保护提供了种群遗传学理论依据，并可能为人类代谢、骨骼肌发育和昼夜节律相关疾病研究提供参考。
Theoretical Framework
N/A
Research Design
本研究属于观察性比较基因组学研究结合功能实验验证。分析单元为基因、基因家族和物种群体。主要自变量为物种（慢懒猴与其他哺乳动物）及基因型（慢懒猴型 vs. 人类型变异），因变量包括：基因家族拷贝数、基因表达水平、酶活性、肌纤维类型基因表达量、蛋白结合能力、有效种群大小等。具体操作化方式：基因家族扩张通过比较基因组学鉴定；正选择通过PAML分支位点模型检测；趋同进化通过氨基酸替换模式比较；功能验证通过体外酶动力学实验、小鼠成肌细胞转染实验和免疫共沉淀实验完成。
Data & Sample
- 基因组组装样本：1只成年雌性倭蜂猴的骨骼肌样本（来自中国科学院昆明动物研究所野生动物种质资源库动物分部），用于PacBio长读长测序、Hi-C测序和转录组测序（8个组织：心、肝、脾、胃、肾、脑、肌肉、舌）。
- 重测序样本：50只倭蜂猴和6只孟加拉蜂猴的全基因组重测序数据。
- 比较基因组学参考：人类、小鼠、猕猴、鼠狐猴、婴猴、考拉、二趾树懒等其他哺乳动物基因组。
- 地理范围：东南亚地区，倭蜂猴与孟加拉蜂猴分布区重叠。
- 时间：样本采集时间未明确说明；种群历史分析覆盖约100万年时间尺度。
Analytical Strategy
- 基因组组装与注释：Canu和WTDBG组装、Quickmerge合并、Pilon校正、Hi-C辅助染色体组装；BUSCO和CEGMA评估完整性；RepeatMasker注释重复序列；多种方法预测蛋白编码基因。
- 比较基因组学：基因家族扩张/收缩分析（CAFÉ类方法）；PAML codeml分支位点模型检测正选择基因；趋同进化分析（检测倭蜂猴与考拉、倭蜂猴与二趾树懒之间的平行/趋同氨基酸替换）。
- 转录组分析：比较不同组织间基因表达，鉴定肝脏特异性高表达基因。
- 功能实验：PITRM1酶活性动力学测定；MYOF转染小鼠成肌细胞后检测快肌纤维基因（MYH2、MYH1）表达；PER2与CRY的免疫共沉淀结合实验。
- 群体基因组学：SNP鉴定（36.3百万个高质量SNP）；系统发育树、PCA、ADMIXTURE群体结构分析；PSMC和MSMC模型重建有效种群大小历史及分化时间。
- 稳健性检验：未明确提及专门的稳健性检验，但通过多方法交叉验证（如BUSCO、CEGMA、短读长比对率）评估组装质量，并通过多物种比较和功能实验验证基因组分析结果。
Results & Findings
基因组组装质量：倭蜂猴基因组组装为2.78 Gb，分布于25条染色体（2n=50），scaffold N50达129.9 Mb，显著优于已有的鼠狐猴和婴猴基因组，为后续分析提供了高质量参考。
解毒能力增强：慢懒猴中多个与杂环化合物结合和代谢相关的基因家族特异性扩张。GSTA基因家族在倭蜂猴单倍体基因组中有14个拷贝，远多于人类、猕猴和小鼠；该家族在肝脏和胃中表达显著高于人类和小鼠。此外，FMO5（I相解毒酶）在倭蜂猴谱系中受正选择且在肝脏中高表达。这些结果支持慢懒猴通过基因扩张和正选择增强了饮食中次生代谢物的解毒能力，与其以树胶/渗出物为主的低热量、高毒素饮食相适应。值得注意的是，婴猴的GSTA家族也显著扩张，提示趋同适应。
PITRM1趋同进化与低代谢率：倭蜂猴与考拉在PITRM1基因上共享3个趋同氨基酸替换（K633R、T697S、D925E），其中两个替换也存在于巽他懒猴中。功能实验表明，倭蜂猴PITRM1的酶活性比人类低约两倍，与小鼠中PITRM1敲低导致基础代谢和耗氧量下降的报道一致，提示该基因功能减弱可能贡献了慢懒猴的低代谢率表型。
肌肉发育相关基因的适应性进化：MYOF、CHRNA1和HMCN2三个与肌纤维发育相关的基因在倭蜂猴中受正选择。功能实验显示，倭蜂猴MYOF上调快肌纤维基因（MyHC2a、MyHC2x）的能力弱于人类MYOF，与慢懒猴慢肌纤维占优的表型一致。此外，倭蜂猴与二趾树懒在NEB基因（编码肌节细丝肌动蛋白结合蛋白）上存在7个趋同氨基酸替换，该基因在骨骼肌中高表达，可能也与缓慢运动相关。
PER2正选择与异常昼夜节律：PER2、ENOX2、PTGER3、PROX1和SIX3五个昼夜节律调控基因在倭蜂猴中受正选择。PER2的S1271A突变被预测影响PER2-CRY结合域，功能实验证实倭蜂猴PER2与CRY1的结合能力显著弱于人类PER2，可能与慢懒猴独特的昼夜节律和冬眠行为相关。
种群历史与物种分化：倭蜂猴与孟加拉蜂猴在约100万年前分化，此后两者有效种群大小呈相反变化趋势，暗示同域分布下的强种间竞争。过去6万年间两个物种的种群规模均大幅下降，可能与人类迁徙或末次冰期寒冷气候有关。基因组分化在染色体端粒区域最为显著，提示端粒区快速进化可能驱动了同域物种形成。
意外发现：GSTA基因家族在婴猴中也显著扩张，提示不同谱系可能独立进化出类似的解毒适应策略；此外，慢懒猴对自身毒素的抗性可能与其增强的解毒能力有关，这一推测虽未直接验证但具有启发性。
Limitations
- 倭蜂猴基因组中注释的蛋白编码基因数量（26,200个）高于其他灵长类，作者承认可能包含因ab initio预测产生的假阳性。
- 功能实验主要基于体外或异源细胞系统（如小鼠成肌细胞），与体内真实生理环境的关联性需进一步验证。
- 趋同进化分析仅基于序列相似性，缺乏对NEB等基因的功能验证。
- 群体基因组学样本量有限（50只倭蜂猴、6只孟加拉蜂猴），可能影响种群历史推断的精度。
- 关于慢懒猴对自身毒素的抗性机制，作者仅提出推测性解释，未进行直接实验验证。
- 样本采集来源单一（主要为圈养或库存样本），可能无法完全代表野生种群的遗传多样性。
Key Contributions
- 首次组装了倭蜂猴的染色体水平基因组（2.78 Gb，scaffold N50 129.9 Mb），显著提升了链猴亚目的基因组资源质量。
- 揭示了慢懒猴解毒相关基因家族（尤其是GSTA）的特异性扩张，为理解其特殊饮食适应提供了分子证据。
- 通过趋同进化分析和功能实验，证明PITRM1酶活性降低与慢懒猴和考拉的低代谢率相关。
- 通过正选择分析和功能实验，证明MYOF功能减弱与慢懒猴慢肌纤维占优的表型一致，NEB趋同进化可能与缓慢运动相关。
- 证明PER2的S1271A突变削弱其与CRY的结合能力，为慢懒猴异常昼夜节律和冬眠行为提供了遗传解释。
- 通过群体基因组学重建了两个慢懒猴物种的种群历史，揭示了同域物种形成后的竞争动态和近期种群衰退，为保护管理提供了理论依据。
- 为人类代谢紊乱、骨骼肌发育异常和昼夜节律相关疾病的研究提供了潜在的候选基因和进化视角。
Key Claims
"We find that many gene families involved in detoxification have been specifically expanded in the slow lorises, including the GSTA gene family with many newly derived copies functioning specifically in the liver." (Abstract)
"Significant decreases in PITRM1 enzymatic activity in these two species may have contributed to their characteristic low rate of metabolism." (Abstract)
"Functional assays demonstrate the decreased ability of one positively selected gene, MYOF, to up-regulate the fast-type muscle fiber, consistent with the lower proportion of fast-twitch muscle fibers in the slow loris." (Abstract)
"The protein product of another positively selected gene in the slow loris, PER2, exhibits weaker binding to the key circadian core protein CRY, a finding which may be related to this species’ unusual circadian rhythm." (Abstract)
"Finally, population genomics analysis reveals that these two extant slow loris species, which co-exist in the same habitat after sympatric speciation, have exhibited an inverse relationship in terms of their demography over the past 1 million years, implying strong inter-species competition after speciation." (Abstract)